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THE MIXING EFFECT OF SUIIFACE WAVES 
Br Chas. S. Slichtek 

Thk following results are based upon the assumption that the "mixing 
effect" when a portion of a liquid is changed in shape is proportional to the 
magnitude of the shear. If we confine our attention to motion in two dimen- 
sions only, the problem before us is to determine to what extent particles of 
a liquid migrate and take up positions adjacent to particles of liquid not 
formerly their neighbors when a prism of liquid of square cross section changes 
to a prism of rectangular cross section. Theory requires that the particles of 
a "perfect" liquid hold the same relative position to each other after change as 
before. The extent to which this is not true in the case of an actual liquid 
constitutes what is designated as "mixing." Thus if the temperature of a cube 
of liquid varies uniformly from one of the faces to the opposite face, then 
after shear the relative distribution of temperature should be precisely the same 
as before shear. The liquid is "mixed" to the extent that the relative distrib- 
ution of temperature has been disturbed. 

If the dimensions of a rectangular parallelepiped of a liquid change from 
1, 1, 1, to 1 + a, 1, 1/(1 + a), the liquid is said to have experienced a 
finite shear of amount a. Hence, according to this definition, unit finite shear 
corresponds to change of cube 1, 1, 1, to parallelopiped 2, 1, 1/2. For rates 
of movement sufficiently low so that turbulent motion does not result, the 
initial assumption made above that mixing is proportional to shear seems suffi- 
ciently obvious if the present definition* of the magnitude of finite shear be 
adopted. 

It is obvious that the amount of mixing that takes place in an element of 
liquid when subjected to unit shear must depend upon properties of the liquid 
which can be readily determined only by experiment. The viscosity, tlie 
density, and various molecular properties probably all play important parts in 
determining the amount of mixing that takes place when an element of a 
liquid changes shape. 

* It should be noted that this is not the definition adopted In the usual treatment of the 
subject of finite shear. 

(170) 
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The primary cause of mixing in liquids, as ordinarily observed, is due to 
the turbulent motion that is set up when the differential velocity from layer 
to layer, or between the liquid and the boundary reaches a certain critical 
amount. It can readily be seen by consulting the diagrams in the papers by 
Hele-Shaw in Nature, Vol. 50, page 34, and Vol. 60, page 446, that slow 
movements in which turbulent motions have been eliminated result in a very 
small amount of mixing, although the amount of shear be relatively large. 
In the wave motion of a liquid on whose surface simple progressive waves are 
maintained, the deformation or strain of any element is periodic and the 
effect of mixing must be cumulative even though the motion be so slow that 
turbulent motion is not set up. The mixing effect of the waves can be ex- 
pressed as a function of the various constants that particularize the special wave 
motion under consideration. It is obvious that the variation of the mixing 
effect with the depth of the particle of liquid below the surface is the first 
problem to be considered. 

The following notation is used in what follows : 

Height of trochoidal wave H in meters 

Wave length X «< " 

Depth of water d '• " 

Distance from the bottom to particle under consideration . . y ♦' '♦ 
The well known theory of trochoidal surface waves enables us to write down 
the following results : 

The radius of the circular orbit of a particle at any depth is : 

r = mH/2 = ^He-^^^v-'^^^, (1) 

in which wi is a proper fraction, the ratio of the radius of the orbit at level y, 
to the radius of the orbit of a surface particle. 

The period of the wave motion in seconds is, if X be measured in meters : 

T=v'2^^=|xi; d^\, (2) 

the formula being correct to about the sixth decimal place. The error remains 
less than 5 per cent for all values of d satisfying* 

Ad > \ <d 

* The approximation consists in omitting tanh iitdfK from the denominator of the theo- 
retical expression for T*. 
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Consider two particles of water which lie in the same vertical line when the 
crest and also when the trough of the wave passes immediately above them. 
Let the radius of the circular orbit of the higher particle be r^ and of the 
lower particle be rj and let the distance between the centers of the two circular 
orbits be a, which may be taken as small as we please. Then the distances 
between these particles are, 

at crest, a + r^ — Vi; 

at trough, a — I'i + Vi. 

Therefore the small rectangular parallelopiped at the crest, of edges 
1 + (^2 - ri)/a, 1, 1/[1 + (rj - r-i)/a] , 

is first sheared to the cube 

1, 1, 1, 

at an intermediate stage, and then this cube is sheared at the trough to the rec- 
tangular parallelopiped 

1 - (^2 - '•!)/«. 1, 1/[1 + {r^ - »-i)/«] ■ 

The total finite shear from crest to trough is therefore of magnitude : 

2(r-2-»-i)/a, 
or, in the limit, 

2dr/dy. 

Therefore the "mixing eiFect" per second of time, or the coefBcient of 
mixing, is 



4 
5 



i; = iAdr/di/ -^ - Xi = 15.6 AH \-le'-(2'-<')A, (3) 



in which A is the mixing effect, or specific coefficient of mixing, due to unit 
shear. This factor must be experimentally determined for each liquid. 

Assume waves of length 10 meters and of height one meter in water 
whose undisturbed depth is 10 meters, then the.rate of mixing at level y is 
given by 

£'=^0 6'*''^^"^'''''". (4) 

in which Hq is the rate of mixing at the surface. 



1911) 



THE MIXING EFFECT OF SUEFACE WAVES 



173 



Table I shows the coefficient of mixing, or the mixing effect for one 
second of time, of a trochoidal progressive wave one meter high and of various 
lengths. The results are expressed in terms of the mixing efiect due to unit 
finite shear. 

TABLE I 



Depth in meters 


Length of waves in meters. 


below surface 


5 


10 


15 


20 


25 


30 


40 





1.310 


0.495 


0.296 


0.175 


0.125 


0.125 


0.062 


1 


0.406 


0.264 


0.176 


0.127 


0.099 


0.077 


0.053 


2 


0.121 


0.138 


0.115 


0.092 


0.077 


0.062 


0.045 


3 


0.032 


0.086 


0.076 


0.068 


0.059 


0.050 


0.040 


4 


0.009 


0.048 


0.048 


0.049 


0.046 


0.040 


0.034 



The above results give the measure of the migration of particles away 
from their initial location or original companions. The numbers may be 
conveniently designated, as above, "coefficients of mixing," although mixing, 
in the common meaning of the term, does not talce place unless another con- 
dition be present, namely, tinless there be a quality varying in magnitude 
throughout the liquid which can be modified by the migrating particles. We 
do not speak of " mixing " a liquid by stirring it unless it has a different quality 
in some parts than in others. Migrating particles can produce no change in 
a perfectly homogeneous liquid. For example, if a liquid be of homogeneous 
temperature throughout, it will remain of homogeneous temperature no matter 
how turbulent the motion. If, however, the initial temperature be not homo- 
geneous (say it falls continuously from the surface downward, as it actually 
does in summer in the water of a lake or pond) then the migration of particles 
discussed above will bring about a new distribution of temperature. The new 
distribution depends as much upon the gradients in the initial distribution as 
it does upon the coefficient of mixing. This problem will be considered only 
in its special bearing upon changes of temperature taking place in a sheet of 
water due to progressive waves upon its surface. We do not aim to consider 
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the mixing due to the turbulent waves near shore, nor to the effects of the 
wind in driving the warm surface water to lee shores. Professor E. A. Birge, 
who has given wide attention to all these matters, fully discusses the subject from 
the physical and practical side in another place.* 

Let there be any quality or characteristic that remains constant in any 
horizontal layer of the liquid but varies in the vertical direction. For ex- 
ample, let there be a vertical distribution of temperature either continuous or 
discontinuous. If the temperature decrease with the depth below the surface 
there will be no tendency to set up convection currents, as long as the tem- 
perature is above 4°. When surface waves are set up in the liquid, the quality 
or characteristic pertaining to the various horizontal layers will gradually 
change, the distribution becoming immediately continuous across the former 
surfaces of discontinuity. Let the initial vertical distribution be defined by 
the equation, 

t =f{y) (5) 

in which y is the depth of any point below the surface of the liquid [the same 
as(y — (Z) above] and in which/represents a single valued function of y. Let 
the coefficient of mixing at any level be <^(.y). At this level the initial gradient 
of < is 

4Hdy^f{y), 
and after mixing it is 

dlldy=\\-^{y)-]f\y). (6) 

Whence the distribution of the quality t after the "first mixing", or first unit 
of time, is given by 

Likewise the distribution of t after the nth interval is given by 

<„ = y[i-</.(y)]»/'(y)^y + c. (7) 

* The results of these extensive studies of the biological, chemical, and physical phenom- 
ena of fresh water lakes is appearing as a series of Monographs published by the Wisconsin 
Geological Survey. The chemical phenomena are treated in Bulletin XXH. The physical 
studies will appear in a bulletin now in preparation. 
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In the important special case in which the mixing is due to waves, we may 

write 

(^(y) = ae-O'J. (8) 

The forms of a and b have been found above. In like manner 

tn = J(l - ae-^)-f'{y)dy + c, (9) 

which permits of simplification for special values ofy. 

Consider the case in which the temperature falls from the surface to the 
bottom as a linear function of y, so that/(y/) = — hy. Then 

tn = -k f(l - a e-^Ydy + c, (10) 

which, for large values of n, becomes 

t„= -k je-"'"'-"' dy = -~\i e-""^", (11) 

in which li stands for the Integral Logarithmus.* 

The constant of integration must be determined by a quadrature. Mechan- 
ical quadrature was used in obtaining the numerical results given below. 

The important case leading to the development of/(y) in a Fourier's 
Series may also be written down. 

Let ' 

/(y) = kj, sin py, 
then, for large values of w 

tn = K^ (- "«)" («'6' +-P')"* (l/«!) sin(i>y - ^). (12) 

«=o 

in which = ta,n-^(ab/p). 

Keturning to the case of equation (11), the Integral Logarithmus per- 
mits numerical results to be obtained as follows: Let 11= 1 meter, \ = 10 



* The Integral Logarithmus is defined by the integral J*J du/logu. See N. Nielsen, 
Theorie des Integrallogarithmus und verwandter Transzendenten, Leipzig, 1900. Since the publi- 
cation (Leipzig, 1909) of Jahnke and 'Emde's Ftmktionentafeln, tables of this and the allied func- 
tions have become readily accessible. Table II was checked by use of the table of Ei.{—x) on 
page 21 of Jahnke and Emde's tables. 
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meters, d = 10 meters, and let n be sufficiently large so that appreciable values 
may be assigned to nA. From previous notation, 

a = 15.6 ^/X^/* ^ 0.484^; 
b = 27r/X = 0.6283 ; e-^ = (0.534)!'. 

As the value of A is not known, the amount of mixing after a given number 
of seconds cannot be found in absolute units but the amount of mixing at given 
relative times can readily be determined, since nA appears as a single variable 
in the results. With the above numerical data, and with nA taken equal to 1/2, 
1, 2, or to elapsed times proportional to 1, 2, 4, the table II has been prepared. 
Table II gives the variation in temperature in a sheet of water 10 meter 
deep, due to the mixing effect of surface waves of wave length 10 meters and 
height 1 meter, the initial temperature being equal to (10 — y), where y is 
the depth below the surface of the water. 

TABLE II 



Elapsed time, n A = 





i 


1 


2 


Value of constant c = 





0.069 


0.13 


0.216 


Depth, y, below 
surface 


Temp, 
degrees 


Temp, 
degrees 


Temp, 
degrees 


Temp, 
degrees 





10 


9.70 


9.45 


9.06 


1 


9 


8.87 


8.75 


8.48 


2 


8 


7.96 


7.92 


7.78 


3 


7 


7.01 


7.01 


6.98 


4 


6 


6.04 


6.07 


6.09 


5 


5 


6.05 


6.11 


5.16 


6 


4 


4.06 


4.11 


4.18 


7 


3 


3.06 


3.12 


3.19 


8 


2 


2.07 


2.12 


2.20 


9 


1 


1.07 


1.12 


1.21 


10 





0.07 


0.13 


0.22 
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Temperature 



The results are best shown graphically as in the figure. Here the original 
temperature is represented by the straight line 0. The modified temperatures 
at time 1/2 , 1 , 2 , 4 , are shown by the curved lines. The diagram is drawn 
for the upper 5 meters only. The curve for the fourth period is sketched only 
approximately, as the computation diiEculties did not warrant greater exact- 
ness. It will be seen that the effect of the waves is gradually to change the 
temperature of the surface layers to a nearly uniform temperature, the abso- 
lute changes in the intervals of time being less than in earlier equal intervals. 
The changes in temperature show a tendency to cease rather abruptly at points 
where the curves cross the original gradient, as at lo, 2q, Sq. These are at depths 
of about 3 meters below the surface for the wave length and the wave height 
under consideration and gradually move downward with the lapse of time. 
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The same diagram may be used to show approximately the mixing effect 
due to waves of twice the height but of same wave length. In this case the 
curve 22 for wave height 2 meters corresponds to the same interval of time as 
the curve 11 for wave height 1 meter. The points 1, 2, 3 move to the left 
and the points Iq, 2^, 3o, move downward at twice the rate in the case of the 
higher waves than in the case of the lower waves. 

Although the mixing effect in the slow motion of a liquid is a magnitude 
of low order, yet there are important cases of motion in which the migrations 
of particles across the lines of flow is very great. An important case occurs 
in the motion of groundwater. It can be shown in this case that the move- 
ment of particles across the lines of flow equals about 30 per cent of the on- 
ward movement in the line of resultant motion. This means, for example, 
that impurities or contaminations are expanded or spread out over wide areas 
during the motion of underground waters. Diagrams showing some experi- 
mental manifestation of this fact will be found in Water Supply Paper No. 140 
of the U. S. Geological Survey, pp. 22-68. The consideration of the mixing, 
or spread of contaminations, in groundwaters lead to the statement of a very 
general problem of much interest, which the writer will consider in another 
place. 

The magnitude of the coefficient A assumed above cannot be even ap- 
proximately assigned with the experimental data now at hand. It seems 
reasonably certain, however, that its magnitude will range higher than even 
the largest of the coefficients of diffusion. 

University of Wisconsin, 
Madison, Wis. 



